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Abstract 

A rigorous method for introducing the variational principle describing relativistic ideal 
hydrodynamic flows with all possible types of discontinuities (including shocks) is presented 
in the framework of an exact Clebsch type representation of the four-velocity field as a 
bilinear combination of the scalar fields. The boundary conditions for these fields on 
the discontinuities are found. We also discuss the local invariants caused by the relabeling 
symmetry of the problem and derive recursion relations linking invariants of different types. 
These invariants are of specific interest for stability problems. In particular, we present a 
set of invariants based on the relativistic generalization of the Ertel invariant. 

Introduction. In this paper we discuss some problems related to ideal relativistic hydro- 
dynamic (RHD) flows in the framework of the special relativity. They are pertinent to the 
description of flows with discontinuities, including shocks, in terms of canonical (Hamiltonian) 
variables based upon the corresponding variational principle and introducing local invariants 
along with recursion relations. These subjects are of interest from a general point of view and 
are very useful in solving nonlinear problems, specifically, nonlinear stability investigation, de- 
scription of the turbulent flows, etc. In particular, the use of the Hamiltonian approach along 
with additional local invariants of the motion and the corresponding Casimirs allows to improve 
the nonlinear stability criteria. The necessity to consider the relativistic flows is motivated by a 
wide area of applications, including the astrophysical and cosmological problems. 



Variational principles for the ideal relativistic hydrodynamic (RHD) flows, without disconti- 
nuities, have been widely discussed in the literature, see, for instance, [1-3] and citations therein. 
As for the nonrelativistic flows, the least action principle is conveniently formulated in terms of 
the subsidiary fields and the corresponding velocity representation known as the Clebsch repre- 
sentation, see [4-9]. These subsidiary fields can be introduced explicitly by means of the Weber 
transformation, [10], see also [3,4]. Alternatively, they naturally arise from the least action 
principle as Lagrange multipliers for necessary constraints. Using these variables allows one to 
describe the dynamics in terms of canonical (Hamiltonian) variables. The nontrivial charac- 
ter of the Hamiltonian approach is due to the fact that the fluid dynamics corresponds to the 
degenerated case, see [11,12]. 

Recently it was shown [13, 14] that the hydrodynamic flows with discontinuities (including 
shocks) can be described in terms of a least action principle, which includes (as well as natural 
boundary conditions) the boundary conditions for the subsidiary fields. In the present paper we 
show that all type of discontinuities can be described by means of the least action principle in 
terms of the canonical variables for the relativistic flows. 



Variational principle. The relativistic least action principle can be formulated in close anal- 
ogy to the nonrelativistic one. We introduce the action A, 

A = J d 4 xC, (1) 

with the Lagrangian density 

C = -e(n,S) + GJ a Q, a , 

G=(l,u B ,e), Q = (</?, /i B ,S), 5 = 1,2,3, 

where vb, 0, <p, fi B represent subsidiary fields; n, S and e(n, S) denote the particle's number, 
entropy and energy proper densities, J a = nu a is the particle current, and u a is the four- velocity, 



u a = m°(1,v/c), u° = 1/yl — v 2 /c 2 ; comma denotes partial derivatives. Small Greek indexes 
run from to 3, and the Latin indexes run from 1 to 3; a; = ct, r ). The metric 

tensor, g al3 , corresponds to the flat space-time in Cartesian coordinates, g a @ = diag{— 1, 1, 1, 1}. 
The four-velocity obeys The normalization condition 

u a u a = g a pu a vP = -1. (3) 

Below we consider the four-velocity and the particle density n as dependent variables expressed 
in terms of the particles current J a , 

u a = J a /\J\, n=\J\ = y/-J a J a . (4) 

The fluid energy obeys the second thermodynamic law 

de = nTdS + n~ 1 wdn = nTdS + Wdn , (5) 

where T is the temperature and w = e + p is the proper enthalpy density, p is the fluid pressure, 
W = w/n. 
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Variation of the action given by Eq. (JTJ) with respect to the variables J a , 6, and Q = 

(ip,fi B ,S), which are supposed to be independent, A = A[J a , <p, fi B , S, v B , 6], results in the 
following set of equations 

5J a : ^Wu a = V a = -GQ, a , (6) 

5<p J a a = 0, (7) 

Sfj, B d a (J a u B ) = 0, or Dv A = 0, (8) 

Su B D/i B = 0, (9) 

<9 a (J Q 6), or L>9 = — T, (10) 

,56 = 0, (11) 

where D = u a d a . Eq. © gives us the Clebsch type velocity representation, cf. Ref. [2]. Con- 
tracting it with u a results in the dynamic equation for the scalar potential ip, 

Dip = W. (12) 

Both triplets fi B and ub represent the advected subsidiary fields and do not enter the internal 
energy. Therefore, it is natural to treat one of them, say, fi B as the flow line label. 

Taking into account that the entropy and particle conservation are incorporated into the set of 
variational equations, it is easy to make sure that the equations of motion for the subsidiary vari- 
ables along with the velocity representation reproduces the relativistic Euler equation. The latter 
corresponds to the orthogonal to the flow lines projection of the fluid stress-energy-momentum 
T a/3 conservation, cf. Ref. [16,17], 

We may then write the relativistic Euler equation as 

(Va,p — Vp^U 13 = TS, a , (14) 



where the thermodynamic relation 



dp = ndW - nTdS (15) 



is taken into account. The vector V a , sometimes called the Taub current, [15], plays an important 
role in relativistic fluid dynamics, especially in the description of circulation and vorticity. Note 
that W can be interpreted as an injection energy (or chemical potential), cf., for instance [16], 
i.e., the energy per particle required to inject a small amount of fluid into a fluid sample, keeping 
the sample volume and the entropy per particle S constant. Therefore, V a is identified with the 
four-momentum per particle of a small amount of fluid to be injected in a larger sample of fluid 
without changing the total fluid volume and the entropy per particle. 
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Boundary conditions. In order to complete the variational approach for the flows with dis- 
continuities, it is necessary to formulate the boundary conditions for the subsidiary variables, 
which do not imply any restrictions on the physically possible discontinuities (the shocks, tan- 
gential and contact discontinuities), are consistent with the corresponding dynamic equations, 
and thus are equivalent to the conventional boundary conditions, i.e., to continuity of the parti- 
cle and energy-momentum fluxes intersecting the discontinuity surface R(x a ) = 0, cf. Ref. [17], 

{J} = 0, J=J a n a , (16) 



{T a %} = 0, (17) 
where n a denotes the unit normal vector to the discontinuity surface, 



N a /N, N a — R >a N= y/N a N* 



and braces denote jump, {X} = X\r =+0 — X|# = _ . 

Our aim is to obtain boundary conditions as natural boundary conditions for the variational 
principle. In the process of deriving the volume equations we have applied integration by parts 
to the term J a G5Q ta . Vanishing of the corresponding surface term along with that resulting 
from the variation of the surface itself will lead to the appropriate boundary conditions after the 
variational principle has been specified. 

Rewriting the (volume) action with the discontinuity surface being taken into account in 
explicit form as 

A = ( d 4 x C?B{sR) , (19) 
J ?=±i 

where 9 stands for the step-function, we obtain the residual part of the (volume) action in the 
form 

&Kes = [d 4 xJ2 [s£5d(R)SR + 6(qR)d a (J a G5Q)} . (20) 

Here 5d denotes Dirac's delta-function and we omit the index q labeling the quantities that cor- 
respond to the fluid regions divided by the interface at R = 0; the superscript q ^ corresponds 
to the quantities in the regions R ^ 0, respectively. Integrating the second term by parts and 
supposing that the surface integral f d A x X] ? =±i ®a {0{qR){u a G5Q)) vanishes due to vanishing of 
the variations SQ at infinity, we arrive at the residual action expressed by the surface integral 

5A \res = [ d ^ X S 5 *>( R ) [ C6R ~ R*J a G5Q] . (21) 

here SQ designates the limit values of the volume variations, 5Q ± = {5Q)r=±o- It is convenient 
to express these variations in terms of variations of the boundary restrictions of the volume 
variables, 5{Xr = ±q) = 5X ± , and variation of the discontinuity surface. It is easy to show that 

5X = 6X+ \N\- l n a X,JR - X^P^bf , (22) 

where = 5 a /3 —n a ni3, and df 13 is an arbitrary infinitesimal four- vector related to the one-to-one 
mapping of the surfaces R = and R + 5R = 0. 



4 



Vanishing of the action variation with respect to variations of the surface variables 5R and 
6fP (which are supposed to be independent) results in the following boundary conditions 

5R :=>• {p+ (u a n a ) 2 w) = 0, (23) 

Sf : => {jyj a iV> 7 } = 0, or P\ { JWu^} = 0, (24) 

which are equivalent to continuity of the momentum and energy fluxes, cf. Eq. (j!7jh Here we 
consider that the 'on shell' value of the volume Lagrangian density, C eq , is equal to the pressure, 
C eq = — e + nGDQ = — e + w = p. 

Now we can complete formulation of the variational principle appropriate both for continuous 
and discontinuous flows. The independent volume variables are indicated above, and independent 
variations of the surface variables are 5R, 5/ 13 , the variations of the surface restrictions of the 
generalized coordinates 5<p, 5fi B , supposed to be equal from both sides of the discontinuity, 
{Sip} = {Sfi B } = 0, and SS with {5S} ^ 0. Under these assumptions we arrive at the following 



subset of the boundary conditions 

5fi { J a n a } = {J} = for {5!p} = 0, (25) 

SJ1 B : => {v B J a n a } = J{v B } = for {6J1 B } = 0, (26) 

6S± J a n Q 9 ± = J6 ± = 0. (27) 



Eqs. (|2~3|) - (|2~H|l reproduce the usual boundary conditions, and Eqs. (|2*7jl are the boundary 

conditions for the subsidiary variables. Other boundary conditions for the latter variables do not 
strictly follow from the variational principle under discussion. But we can find them from the 
corresponding volume equations of motion, providing, for instance, that they are as continuous 
as possible. 1 The natural choice corresponds to continuity of their fluxes, 

{n a u a nfi B } = J{fi B } = 0, (28) 

{n a ii a mp} ee J{^} = 0. (29) 

The set of the boundary conditions given by Eqs. (|23|) -(|29| ) is complete and allows one to 
describe any type of discontinuities, including shocks. For the latter case J ^ and we arrive 
at continuity of the variables Ub, fi B , <fi and zero boundary value of 0. For J = the flow lines 
do not intersect the discontinuity surface and we obtain very weak restrictions on the boundary 

^^Note that the choice of the boundary conditions for the fields ip, fi B , vb and O is not unique due to the 
fact that they play the roles as generalized potentials and therefore possess the corresponding gauge freedom 
relating to the transformations <p, fi B , vb, 6 — > <f', [t' B , v' B , 0' such that u' a = u a (given by the representation 
(jnj). For instance, it seems possible to use entropy S as one of the flow line markers. But if we are dealing 
with discontinuous flows then it is necessary to distinguish the Lagrange markers of the fluid lines. fi B , and the 
entropy, S. Namely, the label of the particle intersecting a shock surface evidently does not change, but the 
entropy does change. Thus, entropy can be chosen as one of the flow line markers only for the flows without 
entropy discontinuities. 
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values of the subsidiary variables, cf. the nonrelativistic case discussed in Refs. [13,14]. Note 
that for the specific case J = (slide and contact discontinuities) we can simplify the variational 
principle assuming all both-side variations of the subsidiary variables to be independent. 

The above variational principle allows modifications. First, it is possible to exclude con- 
straints, expressing the four- velocity by means of the representation (jHJ). In this case the volume 
Lagrangian density can be chosen to coincide with the fluid pressure, cf. Ref. [2], where the 
continuous flows are discussed in detail. Second, we can include into the action the surface term 
respective for the surface constraints, cf. Refs. [13,14,18,19], where such surface terms are 
discussed for ideal hydrodynamics and magnetohydrodynamics in the nonrelativistic limit. This 
can be done for the cases both with excluded and non excluded volume constraints. 

Canonical variables. Starting from the action in Eq. and Lagrangian density given by 
Eq. (j2J) we can introduce the canonical (Hamiltonian) variables according to the general receipt. 
Let Q represents the canonical coordinates then 

8 A 

P = —- = J°G = (7T^TT flB ,n s ) (30) 

provides the conjugate momenta. Relations (|3U|) cannot be solved for the generalized velocities 
(5,o suggesting that we are dealing with the degenerated (constraint) system, cf. Refs. [3,5,11,12]. 
But the constraints are of the first type. Thus, performing the Legendre transform with respect 
to Q we arrive at the Hamiltonian density 

H = PQ, -p(W 7 S) t (31) 

where we suppose that the four- velocity is given by representation (JBJ). Making use of the defi- 
nition (|30|) and of the time component of the velocity representation, Eq. ©, we can transform 
the first term in Eq. (J31j) as 

PQ fi = J°GQ fi = -tv^Vo = n^V°. (32) 
Taking into account the normalization condition for the Taub current, V a V a = —W 2 , we obtain 



V° = ^W 2 + V a V a . (33) 
Consequently, we arrive at the following Hamiltonian density 



H = H(P, Q, Q, a ; W) = ^W 2 + V a V a ir v - p(W, S). (34) 

In terms of the canonical coordinates and momenta the space components of the velocity are 

n v V a = -PQ, a . (35) 

The canonical equations following from this Hamiltonian reproduce in a 3 + 1 form the dynamical 
equations above for the variables entering the Taub current representation. Variation of the 
action with respect to the chemical potential W results in the identity 

n = , ^ . (36) 

y/1 + V a V a /W 2 



6 



Obviously, this relation is equivalent to Eq. (J33)) . expressing the particle density n in terms of 
the variables entering the Hamiltonian. 

We emphasize that the Hamiltonian given by Eq. (}3"4*|) depends not only on the generalized 
coordinates <p, fi B , S, their spatial derivatives and conjugate momenta, but also on the chemical 
potential W as well. Evidently, we can consider W as the additional generalized coordinate with 
zero conjugate momentum, irw = 0. This condition is consistent with the dynamic equations 
due to the fact that d ir w = dH/dW = 0, cf. Eq. (jSfij). 

Bearing in mind that we are dealing with flows having discontinuities, it is seen that in 
the discussed variant of the least action principle we do not arrive at the additional surface 
variables except for the variable defining the discontinuity surface, R. But this enters the 
action functional without derivatives. Therefore, the corresponding conjugate momentum is 
zero-valued. Introducing the Hamiltonian variables for the flows with discontinuities we have to 
treat R as the surface function, defining some (surface) constraint. The latter is nothing else 
than continuity of the normal component of the fluid momentum flux, Eq. (|23|). 

Local invariants and recursion relations. In addition to energy, momentum, and angular 
momentum conservation, for the ideal hydrodynamic flows there exist specific local conserva- 
tion laws related to the advected and frozen-in fields, and corresponding topological invariants 
(vorticity, helicity, Ertel invariant, etc.), cf. Refs. [3-5,9,20] and citations therein for the non- 
relativistic case. They are caused by the relabeling symmetry, cf. Ref. [9]. Discussion of these 
problems along with the recursion relations linking the four different types of invariants for the 
relativistic flows seems insufficient or absent in the literature , see Refs. [1-3, 15, 21] and cita- 
tions therein. Exploitation of the above description permit us considering these invariants in a 
simplified form. Here we shall briefly discuss the invariants and recursion relations. 
The local Lagrangian invariants, say /, correspond to advected (dragged) quantities, 

DI = 0. (37) 

The partial derivative of each scalar Lagrange invariant gives us the simplest example of the 
Lamb type momentum, L a , which satisfy the following relations 

DL a + u^Lp = 0, (38) 

The next type of invariants are vector conserved quantities, X a , being proportional to the four- 
velocity, i.e., 

i; = 0, X a = \X\u a . (39) 

The trivial example of such quantities is the particle current J a = nu a . The last type corresponds 
to the frozen-in fields, M a p, defined as the antisymmetric tensors obeying the following equation 

DM a/3 + u\M^ + tfpMcr, = 0, M Pa = -M aP . (40) 

Now we can derive some recursion relations. First, an arbitrary function of the Lagrangian 
invariants is also Lagrangian invariant, 

r = F(/,/',...). (4i) 
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Second, multiplication of any invariant by the Lagrange invariant results in the invariant of the 
same type. Symbolically, 



Third, 



L' a = IL a , X ,a = IX a , M' ap = IM af3 . (42) 

L a = I, a , (43) 

M a p = L a ,j3 ~ L Pi „ , (44) 

M a(3 = L a L'p - LpL' a , (45) 

L a = M a/3U P , (46) 

X a = e a ^l(]I'J" v , (47) 

/'" = n - l e a ^uJ^rj" v , DI'" = 0. (48) 

Here is Levi-Civita tensor. Note that the latter relation follows from the fact that the 

conserved current defined by Eq. (|47|) is collinear to the four- velocity, X a = \X\u a = — \X\J a /n. 
Consequently, |X| = — X a u a , and \X\/n = I is an advected scalar. Note also that we can easily 
arrive at the conserved four-currents, say Z a , which are not necessarily collinear to the velocity 
field. In analogy with relation (JTTj) we find 

Z a a = for Z a = ^L p J tV . (49) 

Other examples of such conserved currents are as follows 

Z a = e a ^{IL^)^ (50) 

Z« = e a ^{lL^) v , (51) 
Z a = (*M a/3 /) , (52) 

where *M af3 is dual to M a p. 

These recursion relations represent a strict analog of the nonrelativistic ones, cf. Refs. [3,5,20] 
and can be proved by direct calculations. Existence of these four types of invariants is related to 
the dimensionality of the space-time as becomes evident if we remember that they correspond 
to the differential forms of 0th, 1st, 2nd and 3rd order, cf., for instance, Ref. [16]. The recursion 
procedure allows one to study the structure of the invariants. 

It is noteworthy that not all invariants could be obtained by means of the recursion procedure. 
For flows of general type there also exists the Ertel invariant (or the potential vorticity), cf. 
Ref. [21]. The corresponding conserved four-current is of the form 

£ a = ~e a ^u Ptl S,„ = -*uj au S, u , (53) 
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where oup^ is the (Khalatnikov) vorticity tensor, oup^ = V^p — Vp >fl , and %o ai/ is its dual. The 
vorticity tensor obeys the following equation 

Du au + LO a pu p >v + cop v u p >a = T„S >a - T a S :U = n~ 2 {n^j) )U - n t „p ta ) . (54) 

It can be proved that 8 a is divergence-free (the easiest way is to use the above velocity repre- 
sentation) and 

S a = EJ Q , DE = 0, (55) 

where E is a direct generalization of the well-known nonrelativistic potential vorticity. If the 
vorticity tensor would be the frozen-in quantity then conservation of the vector given by Eq. (|53*J) 
would follow from the above recursion relations. But for the non-barotropic flows it is not so. 
Here u av obeys Eq. (JH%|) and becomes a frozen-in field only for the barotropic (isentropic, in 
particular) flows. Nevertheless, the Ertel current is conserved. Therefore, the basic invariants 
may be obtained by direct calculations. 

It is interesting to note that for the non-barotropic flows there exists a conserved current 
generalizing the helicity current. Consider the reduced Taub vector, 

V a = V a + QS >a , (56) 

where G obeys Eq. (|1U|). and the corresponding reduced vorticity tensor 

= V/3,a ~ V a> p = V A,^% ~ "A,aA*$ ■ ( 57 ) 

This tensor is orthogonal to the flow lines, 

u a pu p = 0, (58) 

therefore, the reduced helicity current 

Z a = *u au V v (59) 

is conserved for arbitrary flows, 

<a = u V U)a = -u) u au = -e ^ J ujp 7 u av = 0. (60) 

Here the "thermassy" field can be chosen in such a way that its initial value is zero and thus 
the initial value of the generalized helicity coincides with the conventional one. 

To derive a set of local invariants we have to start with those following directly from the 
hydrodynamic equations and then apply the above recursion relations. As the simplest example 
consider the non-barotropic flow. Then one can start with the specific entropy S and the Ertel 
invariant E. It is easy to show that the general form of the gauge-independent scalar invariants 
is 

I = F(S,E), (61) 

where F is arbitrary function, cf. the nonrelativistic case, [3]. The structure of the complete 
set of invariants differs for the different type of flows and will be discussed in forthcoming 
publications. For instance, for the barotropic flows we have the independent frozen in field 
M a p = u a p in addition to the scalar invariants of the first generation, S and E. This fact allows 
one to obtain a more complicated set of the invariants. 
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Conclusion. We have shown that it is possible to describe the relativistic ideal fluids with all 
physically allowable discontinuities in terms of the least action principle both in the Lagrangian 
and Hamiltonian description. The boundary conditions for the subsidiary variables, entering the 
Clebsch type velocity representation, are obtained in two different ways: one way follows from 
the variational principle as natural boundary conditions while the other one was obtained from 
the dynamical equations under the assumption relating to the absence of the corresponding 
sources and the maximal continuity compatible with the volume equations. It is possible to 
change the variational principle in such a way that all boundary conditions will result from it, 
i.e., they become natural boundary conditions. For this purpose it is necessary to modify the 
variational principle by adding a surface term with corresponding constraints, similarly to the 
nonrelativistic case (compare with the papers [13,14] for the hydrodynamics and [18,19] for the 
magnetohydrodynamics). These variants will be discussed in future works. 

The approach discussed in this paper allowed us to give a simple treatment of the additional 
invariants of motion and present a set of recursion relations linking different types of invari- 
ants. In particular, we presented a generalization of the helicity invariant for the non-barotropic 
relativistic flows. This approach is suitable for the general relativity and for the relativistic mag- 
netohydrodynamics as well. The discontinuous flows for the general relativity can be described 
in analogy with the above discussion and the results will be published elsewhere. 
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